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Abstract 

Using a hadron and string cascade model, JPCIAE, the energy and centrality dependences of 
charged particle pseudorapidity density in relativistic nuclear collisions were studied. Within the 
framework of this model, both the relativistic p + p experimental data and the PHOBOS and 
PHENIX Au + Au data at y^^„=130 GeV could be reproduced fairly well without retuning the 
model parameters. The predictions for full RHIC energy Au + Au collisions and for Pb + Pb 
collisions at the ALICE energy were given. Participant nucleon distributions were calculated based 
on different methods. It was found that the number of participant nucleons, < Npart >, is not 
a well defined variable both experimentally and theoretically. Therefore, it is inappropriate to 
use charged particle pseudorapidity density per participant pair as a function of < Npart > for 
distinguishing various theoretical models. 
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The main focus of the Relativistic Heavy-Ion Colhder (RHIC) at Brookhaven National 
Laboratory (BNL) is to explore the phase transition related to the quark deconfinement 
and the chiral symmetry restoration. The first available experimental data were the energy 
dependence of charged particle pseudorapidity density in Au + Au collisions at v^„„=56 
and 130 GeV from the PHOBOS collaboration After that, the PHENIX collaboration 
published their data of centrality dependence of the charged particle pseudorapidity density 
in Au + Au collisions at V^„„=130 GeV p. 




FIG. 1: The energy dependence of the charged particle pseudorapidity density at mid- 
pseudorapidity in relativistic p + p and A + A collisions. 

It has been predicted that the rare high charged multiplicity could indicate the onset of 
the Quark-Gluon-Plasma (QGP) phase, since the extra entropy in the QGP phase could 
manifest itself as a huge number of produced particles in the final state 0, ||. On the 
other hand, in |^ the centrality dependence of charged multiplicity has been proposed to 
provide information on the relative importance of soft versus hard processes in particle 
production and therefore provide a means of distinguishing various theoretical models for 
particle production. 

The pQCD calculation with assumption of gluon saturation |^ (referred to as EKRT 
model later) was first used to study the centrality dependence of the charged particle pseu- 
dorapidity density at RHIC. In |^ the HIJING model met with success in describing both 
the energy and centrality dependence of the charged particle pseudorapidity density. The 
conventional eikonal approach and the high density QCD (referred to as KN model later) 
|§| were also used to investigate the centrality dependence and the both methods surpris- 
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ingly obtained almost identical centrality dependence. Recently, authors in reported 
their results from the Dual Parton Model. It was found that the experimental observation, 
the charged particle pseudorapidity density per participant pair slightly increasing with 
< Npart >, was reproduced by [^, |, y], but contradicted the results of 

In this letter a hadron and string cascade model, JPCIAE |]TD[, was employed to study 
this issue further. Within the framework of this model the experimentally measured energy 
dependence of the charged particle mid-pseudorapidity density per participant pair both 
in relativistic p + p and Au + Au collisions at RHIC was reproduced fairly well without 
retuning the model parameters. The predictions for the full RHIC energy Au + Au collisions 
and for Pb + Pb collisions at the ALICE energy were also given. In studying centrality 
dependence the focus was put on the calculations of < Npart >, its definition and uncer- 



tainty. Both the PHENIX [0 and the PHOBOS observations that the charged particle 
mid-pseudorapidity density per participant pair slightly increases with < Npart > could be 
reproduced fairly well by JPCIAE. However, this study indicated that it is not suitable to 
use the charged particle mid-pseudorapidity density per participant pair as a function of 
< Npart > to constrain theoretical models for particle production, because < Npart > is not 
a well defined physical variable both experimentally and theoretically. 

The JPCIAE model was developed based on PYTHIA Q. In the JPCIAE model the 
nucleons in a colliding nucleus are distributed randomly in the sphere of the nucleus with 
a radius of 1.12^4^/^ fm. The modules of the nucleons are sampled by the Woods-Saxon 
distribution and the solid angles of the nucleons are sampled uniformly in An. Each nucleon 
is given a beam momentum in z direction and zero initial momentum in x and y directions. 
After the construction of initial particle list the collision time of each colliding pair is calcu- 
lated under the requirement that the least approaching distance of the colliding pair along 



their Newton straight-line trajectory should be smaller than yatot/T^, where atot refers to 
the total cross section. The nucleon-nucleon collision with the least collision time is then 
selected from the initial collision list to perform the first collision. After the first collision, 
both the particle list and the collision list are updated and now the collision list may consist 
of not only nucleon-nucleon collisions, but also collisions between produced particles and the 
nucleons and between produced particles themselves. The next collision is selected from the 
new collision list. The processes proceed until the collision list is empty. 

For each collision pair, if its CMS energy is larger than a given cut, we assume that 
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strings are formed after the collision and PYTHIA is used to deal with particle production. 
Otherwise, the collision is treated as a two-body collision p3| , |14|, |1^]. The cut (=4 GeV in the 
program) was chosen by observing that JPCIAE correctly reproduces charged multiplicity 
distributions in AA collisions. 

It should be noted here that the JPCIAE model is not a simple superposition of nucleon- 
nucleon collisions since the rescatterings among participant, spectator nucleons and produced 
particles are taken into account. We refer to [l^ for more details of the JPCIAE model. 
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FIG. 2: The number of participant nucleons < N-part >, (a), and the ratios of the different curves 
in (a) to the PHENIX data, (b), as a function of the percentage of geometrical cross section. 



In Fig. 1(a) the experimental data of charged particle pseudorapidity density per partic- 
ipant pair at mid-pseudorapidity in relativistic p + p (open triangles and rhombuses with 
error bar) and A + A collisions (open circles with error bar for Pb + Pb at SPS and full circles 
with error bar for Au + Au at RHIC) [|T| are compared with JPCIAE model (full stars for 
p + p and the solid curve for A + A collisions), HIJING model (dotted-dash curve without 
jet quenching and dashed curve with jet quenching) ^ and EKRT model (dotted curve) 
1^. The data of both p + p and A + A collisions at relativistic energies were reproduced 
fairly well by JPCIAE model without retuning model parameters. Fig. 1 (b) is the same 
as (a) but the vertical coordinate here is the charged particle pseudorapidity density itself. 
The JPCIAE model predictions for full RHIC energy Au + Au collisions and for Pb + Pb 
collisions at the ALICE energy in both panels may supply a benchmark for QGP formation 
since the QGP phase is not included in the JPCIAE model. 

Since number of participant nucleons, < Npart >, plays a crucial role in the presentation 
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of PHOBOS or PHENIX centrality dependence data we first make a study on < Npart >■ 
In the fixed target experiments the participant nucleons from the projectile nucleus with 
atomic number A, for instance, is estimated from 



Npart - A* {1- Iff) 
beam 



(1) 



where Ezdc refers to the energy deposited in the Zero Degree Calorimeter dominated by 



projectile spectator nucleons and E^l^^ is the kinetic energy of beam [[l^]. However, in 
the collider experiments, in order to obtain < Npart > one has to relate the measurables 
to Monte Carlo simulations. In PHENIX, for instance, simulations for the response of the 
Beam-Beam Counter and the ZDC were used to calculate < Npart > via a Glauber model 
0. In PHOBOS < Npart > is derived relating HIJING simulations to the signals in the 
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FIG. 3: The charged particle pseudorapidity density at mid-pseudorapidity in Au + Au collisions 
at -y/s„„=130 GeV as a function of the percentage of geometrical cross section. 

On the theory side, first, < Npart > could be calculated geometrically (referred to as 
method a later) |jT^, number of participant nucleons from the projectile nucleus, for instance, 
reads 

N^artib) =poJ dVe{Rp - [x^ + {h- yf + z^\^l^)d{RT - {x^ + y")"^). (2) 
Second, in the Glauber model Npart is calculated by (referred to as method b later) 0] 
Npartib) = / d^sTA{b - s)[l - exp{-ainTB{s))]+ / d^sTB{s)[l - exp{-ainTA{b - s))], (3) 



where am ~40 mb is the inelastic nn cross section at RHIC and Ta refers to the nuclear 
thickness function of nucleus A. The third method is to count the participant or the spectator 
nucleons in the simulation for nuclear collisions and then to average over simulated events. 
However, there is multifarious in simulating programs and even in the definition of the 
participants and spectators. FRITIOF 7.02 (referred to as method d later) |T^ was popularly 
employed in the past. In FRITIOF the wounded nucleons, i.e., nucleons which suffer at least 
one inelastic collision, are counted and identified as < Npart >■ It should be pointed out here 
that in FRITIOF leading nucleons undergo multiple scatterings and get excited (forming 
strings) during collisions, but produced particles from the string fragmentation do not have 
rescattering. Unlike JPCIAE, FRITIOF is not a transport model, there is no space-time 
coordinates associated with each particle. In JPCIAE simulations we have devised three 
counting methods: First, the leading nucleons involved in at least one inelastic nucleon- 
nucleon collision with string excitation are counted and identified as < Npart >■ This is 
called method c. It should be mentioned that the JPCIAE results in Fig. 1 were calculated 
by < Npart > from method c. Second, in method e, spectator nucleons are counted at the 
final state of JPCIAE simulations without rescattering (i.e., only nucleon-nucleon collisions 
are included), then < Npart > is calculated through 

< Npart >={A + B)- < Nspec >, (4) 

where A and B refer to the atomic numbers of the target and projectile nuclei. Third, 
method f is the same as method e but JPCIAE simulations are with rescattering. The 
difference among method c, e, and f is that those nucleons which only experience two-body 
nucleon-nucleon collisions (without string formation) are included into < Npart > in method 
e, while in method f even those nucleons which suffer collisions with other produced particles 
are also included into < Npart >■ In other word, in method e and f the nucleons that are 
knocked out of the colliding nuclei by the produced particles are included into < Npart >■ 
In emulsion chamber experiment, such nucleons (protons) are usually called 'grey tracks'. 

Fig. 2 (a) gives < Npart > calculated by different methods as a function of percentage 
of the geometrical (total) cross section and compares them with results of PHENIX (solid 
circles with error bar) [0 and of PHOBOS (open circles with error bar) |Tl|. The thin 
dotted, dashed, and sohd curves in this panel are the results of method a (geometry), 
b (Glauber model, taken from |^), and c, respectively. Thick dotted, dashed, and solid 
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curves are, respectively, the results of method d (FRITIOF), e, and f. One knows from 
Fig. 2 (a) that except method f < Npart > from different methods are close to each other 
(the 10%— 15% difference should contribute to the systematic error of the experimentally- 
extracted < Npart > ) ioT most central collisions but the discrepancies among them increase 
with decreasing centrality in general. It is surprising that the results of geometry method 
are closest to the results of PHINEX or PHOBOS. In Fig. 2 (b) the ratios of < Npart > 
from methods a, b, c, d, e and f to the corresponding results of PHENIX are given. 

The charged particle pseudorapidity density at mid-pseudorapidity in Au + Au collisions 
at A/i„„=130 GeV as a function of the percentage of geometrical cross section is given in 
Fig. 3. In this figure the full and open circles with error bar are the PHENIX and 



PHOBOS |[TT[] data, respectively. The full and open triangles, respectively, are the results of 
JPCIAE model with rescattering and without rescattering. The abscissa of the data point 
in this figure is set at the middle of the corresponding bin of percentage of the geometrical 
cross section. Globally speaking, the experimental data were reproduced fairly well by the 
JPCIAE model. However, the agreement between the experimental data and the JPCIAE 
model becomes less satisfied for peripheral collisions. One knows from this figure that the 
rescattering only leads to a few percent increase in charged multiplicity although rescattering 
might enhance yields of strangeness, S~ + S~ for instance, by a couple of times. 
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FIG. 4: The charged particle pseudorapidity density per participant pair at mid-pseudorapidity 
in Au -|- Au collisions at y^„^=130 GeV as a function of the number of participant nucleons, 

< Npart >■ 



In panel (a) of Fig. 4 we compare the PHENIX data of the charged particle mid- 
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pseudorapidity density per participant pair (full circles with shaded area of systematic errors) 
with JPCIAE model (full triangles, < Npart > from method c) and with HIJING model 
(dotted curve), KN model (solid curve), and EKRT model (dashed curve). One can see that 
except EKRT, all the other three models predict an increase of {drich/ dri\^=Q) / (0.5 < Npart >) 
as a function of < Npart > though the theoretical results seem to underestimate the PHENIX 
data. Fig. 4 (b) compares the PHENIX data to the JPCIAE results of < Npart > calculated 
by method a (thin dotted curve), b (thin dashed curve), c (thin solid curve), d (thick dotted 
curve), and e (thick dashed curve), respectively. One sees from this panel that starting from a 
single result of charged particle mid-pseudorapidity density from the JPCIAE model but us- 
ing < Npart > from different methods, it is possible to lead {dnch/dri\rf=o)/ (0.5 < Npart >) to 
either increase or decrease with the increase of < Npart >■ Although < Npart > from method 
a are closest to the PHENIX results (cf. Fig. 2 (b)) the {dnch/dri\r^=o) / {0.5 < Npart >) from 
JPCIAE have actually centrality dependence opposite to the PHENIX result (cf. Fig.4 (b)) 
because the {drich/ dri\r^=Q) from JPCIAE is higher than the PHENIX result for peripheral 
collisions (cf. Fig. 3). On the other hand, even though the discrepancy between < Npart > 
from method c and PHENIX slightly increases with decrease of centrality (cf. Fig. 2 (b)) 
the {dnch/dr]\n=o)/{0.5 < Npart >) from JPCIAE with method c is close to the PHENIX 
data. If dnch/dT]\rj=o from EKRT model was not normalized by < Npart > from method b, 
as did in [0, but by < Npart > from method d the results of {drich/ dri\rf=o) / {0.5 < Npart >) 
might have somewhat similar centrality dependence as in PHENIX data. Therefore one 
learns here that it is hard using {dnch/dri\rf=o)/ (0.5 < Npart >) as a function of < Npart > to 
distinguish various theoretical models for particle production since < Npart > is not a well 
defined physical variable. 

In summary, we used a hadron and string cascade model, JPCIAE, to investigate 
the energy and centrahty dependences of charged particle pseudorapidity density at mid- 
pseudorapidity in relativistic p + p and A + A collisions. Both the relativistic p + p exper- 
imental data and the PHOBOS and PHENIX data of Au + Au collisions at RHIC could 
be reproduced fairly well within the framework of the JPCIAE model without retuning any 
parameter. The JPCIAE model predictions for full RHIC energy Au + Au collisions and for 
Pb+Pb collisions at the ALICE energy are also given. This study shows that since < Npart > 
is not a well defined physical variable both experimentally and theoretically it may be hard 
to use charged particle pseudorapidity density per participant pair at mid-pseudorapidity as 
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a function of < Npart > to distinguish various theoretical models for particle production. 
Finally, the financial supports from NSFC in China and DOE in USA are acknowledged. 
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